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We study the effect of temperature shift on aging phenomena in the Random Energy Model 
(REM). From calculation on the correlation function and simulation on the Zero-Field-Cooled 
magnetization, we find that the REM satisfies a scaling relation even if temperature is shifted. 
Furthermore, this scaling property naturally leads to results obtained in experiment and the 
droplet theory. 
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§1. Introduction 

The aging, which are dynamical behaviors largely depending on the history of system, is one 
of the most striking phenomena in the complex systems such as spin glasses, glasses, polymers 
and proteins, and has been studied vigorously from both theoretical and experimental aspects. In 
experiments on spin glasses, one of the most familiar methods for the study of the aging phenomena 
is the m ea_, of the ZeroField-Coo.ed (ZFC) — atio Jft In this _ent, the 
sample is quenched to a temperature T below T c in zero field. After a waiting time i w , a weak 
magnetic field is applied and the magnetization M{t) is recorded as a function of observation time 
t. In this experiment, the length of t w determines the degree of equilibration and aging appears as 
the t w dependence, especially, the peak of the relaxation rate S(t) = emerges around t w as if 

the system remembers how the equilibration has proceeded before the observation. 

Granberg et al$ modified this experiment by shifting the waiting temperature from T to T — AT 
and observed that, although t w was unchanged, the maximum of S(t) shifted to left with increasing 
AT as if t w decreased. We refer to the position of this maximum as apparent wait time t^ p . This 
suggests that the equilibration at T — AT during t w only corresponds to the one at T during t^ pp . 
For AT < 0, the position shifts to right. They also investigated the relation between t w and AT 
and found 

log(C P Aw) ~ ~a(t w )AT, (1.1) 
with the result that the coefficient a(t w ) is a monotonically increasing function. 
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From the theoretical aspects, this phenomenon is demonstrated by the droplet theory very well,tt 
and is reproduced in the study of the two and three dimensional EA ising spin glass mo del& and 
a hierarchical diffusion model® In this manuscript, we study this phenomenon with the Random 
Energy Model (REM)JJiB At first, we calculate a correlation function with the temperature 
shift and confirm that this correlation function satisfies a scaling law. Next, we carry out the 
simulation to observe the temperature shift ZFC magnetization. As a result, we confirm that this 
magnetization also satisfies the same scaling law as the correlation function, which naturally leads 



to the relation eq. (1.1). 

The organization of this manuscript is following. In section 2, we introduce the REM. In sec- 
tion 3, we define the temperature shift correlation function and calculate it. In section 4, we 
observe the temperature shift ZFC magnetization by simulation and analyze the mechanism of this 
phenomenon. In section 6, we summarize this paper with some discussions. 

§2. The Random Energy Model 

The REM is schematically shown in Fig. l.B The bottom points represent the accessible states 
of this system. The length of each branch represents the barrier energy E, over which the system 
goes from one state to another. This energy is an independent random variable distributed as 

d E 

p(E)dE=—exp[-E/T c ), (2.1) 

where T c is the transition temperature. 

From the Arrhenius law, the relaxation time r(a : T), the average time for the system to escape 
from the state a at T, is related to E(a) as 

r(a : T) = tq ex.p[E(a)/T], 

where To is a microscopic time scale. Hereafter To is used as the unit of time and set to 1. From 
eq. fl2.1| ), the distribution of r can be written as 

Px(r)dr = -^dr (t>1), (2.2) 

where x = T/T c . From eq. ( |2.2D , it is easily shown that the averaged relaxation time (t) is x/(x — 1) 
for x > 1 and infinite for x < 1. This means that the transition from the ergodic phase to the 
non-ergodic phase occurs at T c . 

For dynamics, we employ the following simple Markoff process. First, the system is thermally 
activated from a in unit time with the probability 

W(a : T) = ex.p[-E(a)/T]. (2.3) 

Since the time evolution is Markoffian, we can easily evaluate the probability q(a,t : T)dt that for 
the event to occur during t + to and t + to + dt with knowing the system is in the state a at to as 

q(a, t : T)dt = W(a : T) exp[-W(a : T)t]dt. (2.4) 
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We generate an event time t according to q(a, t : T) to perform event-driven Monte Carlo simulation. 
When the magnetic field H is applied, we replace E(a) in eq. ( |2.3j ) with E + HM a , where M a is 
the magnetization of the state a. The magnetization is chosen randomly and independently from a 
distribution D{M) with zero mean. In this simulation, we choose the uniform distribution between 
-1 and +1 for D(M). 

After the activation, the system falls to one of all states with equal probability. In the limit 
N — > oo, we can neglect the possibility that the system has stayed one state more than twice in our 
finite time simulation. Therefore, in this simulation, we create a new state by using the distribution 
functions p(E) and D(M) whenever activation occurs. 

§3. The Temperature Shift Correlation Function 

At first, we define the correlation function mentioned in §|]. This correlation function is observed 
in the following procedure. We consider the infinitely high temperature limit for the initial con- 
dition. In other words, the initial state is chosen randomly at t = 0. Then, the system is kept at 
T — AT during t w and heated (or cooled) to T after that. Here, the magnetic field H is always 
zero. In this procedure, we observe the correlation of magnetization between the time t w and t w + t, 
which we hereafter refer to as Cat(^w). It is explicitly written as 

C AT (t,t w ) =J2 M pG /3a (t)M a P a (t w ), (3.1) 

a,f3 

where P a (t w ) is the probability that the system is found at state a at t w and G/3 a {t) is the probability 
that the system which initially stays at a reaches (3 at t. Because M a is independent random variable 
with zero mean, eq. ( |3.1| ) is reduced to the autocorrelation function by taking the average over the 
distribution of the magnetization: 



Cat(Mw) = M 2 Y,G aa (t)P a (t w ), (3.2) 

a 

where M 2 is the variance of the distribution D{M). Finally, we ignore the possibility that the 
system is activated from a between t w and t + t w and we still find the system at a at time t + t w 
in the limit N — > oo. As the result, eq. ( |3.2| ) is further reduced to 



C AT (t, t w ) = M 2 ^2 exp[-W(a : T)t]P a (t w ). (3.3) 

a 

Bouchaud and Dean!) calculated the Laplace transformation of P a (t w ) an d obtain the result, 

p (F) _ E X1 r(a:T-AT) 
a{ ' Nx 1 c(x l ){Er(a:T-AT) + iy 1 ' 

where xi = and 

c(x) = T(x)T(l - x). (3.5) 
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The inverse Laplace transformation of eq. (3.4) leads us to 



Pa {tv 



dss Xl 1 exp 



iVa:ic(xi)r(xi) Jo 
Now, let us rewrite W(a : T) as 

W{a : T) =t(q : T — AT) -7 



r(a:T- AT) 



where 7 



T—AT 



Substituting this equation and eq. (3.6) into eq. (|3.3|) , we obtain 



Cat(£, t w ) 



M 2 



drr Xl i exp[-r 7 t] / dss Xl 1 exp 



(3.6) 



(3.7) 



(3.8) 



c(xi)r(xi) Ji 

By changing the variables asu = t w /T,v = s/r and replacing the upper limit of u- integral with 00 
in the assumption t w ^> 1, we find the following scaling relation, 



Cat(Mw) = C&r(t/ts), 



where 



and 



Cat(X) 



M 2 



00 du 



exp[— Xu 1 — u] 



u dv 
v 



v xi e v . 



(3-9) 

(3.10) 
(3.11) 



c(xi)r(xi) Jo u 

We can easily check that for the case AT = 0, this relation is reduced to the t/t w scaling as shown 
in ref. [7|. 

Now, we are interested in the asymptotic behaviors of Cat(^) in the limit X < 1 and I>1. 
For X ^> 1, we obtain 



M 2 



1 



ds 



exp[— s — [s/X)~r] 



7c(xi)r(xi) Jo s 

—— - / — e- s (s/X) x ' 

■i)r(*i) Jo « v 7 ; 
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xi7c(xi 

r(x 2 ) 



xi7c(xi)r(xi 



-X" 



(3.12a) 



where X2 = y-. In the case X <C 1 and xi + 7 > 1, we find 



1 



7c(xi)r(xi) Jo s 



50 ^ exp [_( s/X )T]{l-e- s } 



(s/X)7 du 



1 



7c(xi)r(xi) Jo s 



ds 



{i-e- s }( s /xy 



(1 - xi)c(xi)r(xi) 



X' 



(3.12b) 
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For ICl and x\ + 7 < 1, 



M 2 c(a:i)r(xi) Jo s 7 « 



:V ^^l+7 e -«(l- S ) 



c(xi)r(xi) Jo s Jo u 
c(x 1 + 7) 



c(xi) 



X. (3.12c) 



Next, we evaluate eq. (3.3) by simulation to check the validity of these results. Random av- 
erage is taken over 10 6 samples. In this simulation, we fix x\ to 0.4 and change X2 and t w as 
0.2,0.4,0.6,0.8,1.2 and 10, 10 2 , 10 3 , 10 4 , 10 5 , respectively. In Fig. 2, the scaling plot for all x 2 is 
shown for (a) t/t$ S> 1 and (b) t/t$ <C 1. We can see that the scaling holds very well and the expo- 
nent in eqs. ( p.!2a ), ( |3,12b| ) and Q3.12cj ) are rather correct, although a little deviation is observed 
on the coefficients for t/t$ <C 1. 

§4. The Result of the Temperature Shift ZFC Simulation 

In this section, we show the result of simulations on the temperature shift ZFC magnetization 
with T c = 1.0 and T = 0.5. The number of samples for random average is 5 x 10 7 . The am- 
plitude of field is 0.1. We prepare the initial condition in the infinitely high temperature limit, 
the same as the temperature shift correlation function. During t w , the temperature of the sys- 
tem is kept at T — AT. After that, the temperature is changed to T and the field H is applied. 
We performed the simulations for the 36 different cases with t w = 10 3 ,3 x 10 3 ,10 4 ,3 x 10 4 and 
AT = 0.08, 0.06, . . . , 0, -0.02, . . . , -0.08. 

In Fig. 3, we plot the relaxation rate S(t) for t w = 3 x 10 and all AT. We can see the position 
of the peak shifts to left with increasing T — AT, although its shape becomes broader. Next, we 
checked whether the same scaling relation as Cat(^w) holds or not. In Fig. 4, we plot S(t) as 
a function of t/t$ for all i w and AT. We can confirm the scaling law. Although the position of 
peak depends on AT to some extent in our simulation, we consider that this dependence vanishes 
if we take AT/T as small as that used in experiment.^ By neglecting this AT dependence, we can 

T-AT 

roughly estimate the position of peak as i^ pp ~ t$ = t w T and obtain the relation eq. ( |1.1[ ) with 
a(t w ) = log (M . This result is the same as the one from the droplet theory including the coefficient 
a(t w )J> 

In order to investigate how the equilibration proceeds during t w , we examine the energy distribu- 
tion P(E,t)(t < t w ) which is defined as the probability density that the system is found at time t 
in one of the states whose energy is E. In Fig. 5, we show how P(E, t) changes with time. At first 
sight, we can see that P(E,t) has a peak at some point E* , which shifts to right with increasing t. 
The value of E* is roughly determined how far the system can be activated during time interval t, 
and estimated as 

T*«Tlogt. (4.1) 
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For E < E* the distribution is well equilibrated, so that the exponent a\ is given as 

Q ' = j ~ %■ (4 - 2) 

Note that this exponent changes its sign at T = T c . While, the other part E > E* is not equilibrated 
and the exponent 02 is equal to that of p(E), i.e., 

oil = ~ (4.3) 

1 c 

Equation Q4.1| ) suggests that the speed of the peak shift varies with the temperature. In Fig. 6, 
we show P(E, t) for t = 10 4 and all AT. We can see that the position of the peak shifts to right 
just as Fig. 3. From these nature of P(E, t), we understand the relation eq. as follows: During 
t w at T — AT, the peak shifts to E* « (T — AT) log(i w ). We can expect that, in this experiment, 
the magnetization most strongly changes at the time scale in which the system can be activated to 
E* at T. Therefore, we can estimate t^P p from the condition 

(T-AT)log(t w )«Tlog(C P ), (4-4) 
which leads to eq. ( |1.1[) with a(t w ) = lc W w ) ; as mentioned above. 

§5. Summary and Discussions 

We have shown that the REM satisfies a scaling law even if the temperature is shifted and time 
1 at rn 
t is scaled by ig = i w T . For the case AT = 0, i/i w -type scaling is reported in the experiment^ 

and simulation on the three-dimensional ±J Ising spin-glass model0> although some systematic 

deviations from the scaling remain on the experiment. It is interesting to check whether these 

systems satisfy the scaling like the REM or not. 

Next, we consider why the droplet theory leads to the same result as the REM. In the droplet 

theory, the crossover from equilibrium relaxation to non-equilibrium one occurs when the droplet 

created after the field application grows to the size of the droplet created before that, and we can 

expect that peak of S(t) emerges there. Therefore, t^ p is determined by the condition 

L(t w ,T-AT)^L(CV,T), (5.1) 

where L(t, T) denotes the characteristic size of the droplet after the age t at T, and is given from 
the condition 

B(T,L(t,T)) = Tlog(t), (5.2) 

and B(T,L) denotes the barrier energy that the system has to overcome to create the droplet of 
size L, and is given as B(T, L) = V{T)L^ explicitly. If T dependence of V(T) is week enough, the 
condition eq. (||) is reduced to eq. ( |4.4| ). 

Recently it was reported that converges to for somewhat large |AT| (|AT/T| > 0.07 in this 
experiment) regardless of its sign© This suggests that the equilibration at T — AT does affect the 
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one at T. In the droplet theory ,E3t§ it is considered that the chaotic change of the equilibrium 
spin configuration against the temperature variation causes this phenomenon. But the REM has 
no reason why this phenomenon is observed. Actually, this phenomenon is not observed in our 
simulation even in the region \AT/T\ < 0.16. But we can show that this situation is improved by 
treating the Multi-layer Random Energy Model (MREM), the details of which will be presented 
elsewhere .111* 
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FIGURE CAPTIONS 



Fig. 1 Structure of the Random Energy Model. 



Fig. 2 (a) Cat(Mw) and (b) 1 — Cat(Mw) is plotted as a function of t/ts for all t w and x<z- Each 
line is the asymptotic behavior estimated in eqs. djjgg ), ( p^b|) and flOgg ). 



Fig. 3 S(t) vs. t for t w = 3 x 10 4 and AT = 0.08, 0.06, . . . , 0, -0.02, . . . , -0.08 (from left to right). 
Fig. 4 The scaling plot of S(t) for all t w and AT. 

Fig. 5 P(E, t) for T = 0.5 at t = 10 5 , 10, ... , 10 3 - 5 (from left to right). 

Fig. 6 P(E, t) for t = 10 4 and AT = 0.08, 0.06, . . . , 0, -0.02, . . . , -0.08 (from left to right). 
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